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We examine the convergence properties of the 2-nucleon EfTective Range Expansion as used in 
Effective Theories (ET-ERE's) for 3-nucleon calculations. We accomplish this by accounting for the 
2-body dynamics with a simple rank-1 separable 2-body potential where the finite range effects can 
be incorporated systematically in both the 2- and 3-body problems. We make our initial comparisons 
in the simple context of the ^Si 2-nucleon channel and the ''S'3/2 3-nucleon channel. We find that 
convergence problems for some of the 3-nucleon scattering amplitudes using the ET-ERE can be 
traced to its poor account of finite range effects that soften the momentum dependence of the 
deuteron propagator in the Faddeev kernel. In contrast, our simple separable potential with dipole 
form factors works very well in all cases considered. 

PACS numbers: 21.45.-v,25.40.Dn,34.50.-s 

I. INTRODUCTION 

Much attention has recently been paid to the low-energy physics of the neutron-deuteron (nd) system using the 
powerful methods of Effective Theories (ET's) [H 121 [11 IH |5T6] . The basis of these approaches is the assumption that, 
if few body systems display long distance features, these should not depend on details of the short distance physics. 
More specifically, most treatments are based on a field theory specified by a Lagrangian which contains a series of 
contact interactions with coefficients depending on increasing powers of momentum. This series, it is hoped, may be 
truncated at some low order and still account for the low energy physics once the (few) parameters of the ET are 
fixed by external inputs. The power of the approach is that these few inputs yield a theory which can then be used 
to calculate systematically any low energy property of the system. 

ET's are very powerful in the limit that the 2-body scattering length is much larger than any other length scale 
in the problem. It is clear that, in this case, microscopic details of the 2-body physics are irrelevant and that any 
solutions found are "universal" in nature; i.e., are applicable to any system for which the limit is valid. In such cases, 
universality arises in the ET's because only the lowest order contact interaction is needed. It is naturally of interest 
to extend this powerful formulation to problems for which the scattering length is large, but not overwhelmingly so 
compared to other length scales characterizing the 2-body interaction. A common means of identifying the appropriate 
2-body length scales is the Effective Range Expansion (ERE) [7 , namely 

pcotS^^- + '^p^ + Vip\.., (1) 
a 2 

where p is the 2-body relative momentum, S is the scattering phase shift at momentum p, a is the scattering length. 
To is the effective range, 174 is the shape parameter, and higher order terms are implicit. When applying ET's to 
low energy phenomena, it is basically the effective range which governs modifications to universal results as one 
leaves the universal regime. Such modifications may be referred to generically as "finite range effects", whether or 
not they can be wholly accounted for by tq alone. In the present work, we examine finite range effects in the low 
energy 3-body problem, focussing on quartet nd scattering, to which ET's have been applied repeatedly over the last 
decade [HI [9l [TOl [iTJ [12] . Several of these previous works have addressed finite range effects in a manner tied closely 
to the ERE as discussed above. 

Here we adopt a different approach which incorporates the finite range of the 2-body interactions at the outset. More 
specifically, we treat the 2-body interaction as a sum of separable terms [H] , each consisting of the outer product of 
form factors which reflect the structure of interactions which are directly related to real or virtual 2-body states. The 
use of separable potentials leads to considerable simplification. Specifically, separable potentials admit an algebraic 
solution to the two-body problem greatly simplifying the three-body integral equations. For the specific purposes of 
this study, the effective range expansion is seen as the zero-range limit of the separable expression whose expansion 
in the separable potential's range parameter illuminates the role of finite range effects allowing a critical examination 
of the manner in which effective range expansions account for finite range effects in the three-body context. 

Separable potentials have a long history in nuclear physics going back to the period immediately following Faddeev's 
seminal work on the 3-body problem (see Refs. [T4j and [T5| for comprehensive lists of representative papers) when 
calculations of Nd scattering relied on finite-range separable NN interactions to an extent that was as nearly universal 
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as the use of contact interactions in ET's is today (see, for example, Refs. [TSl UHl [IZ] ) • Use of this approach continues 
today, though on a much more limited basis. Examples in nuclear physics include References [T8j [19j [20l [21]. The 
method has also been applied quite recently to the new problem of cold, dilute atomic gases in, e.g., References 

In the present work we exploit separable 2-body interactions to illuminate modern treatments of 3-body systems 
using effective theories which begin by assuming contact interactions and then rely on the Effective Range Expansion 
to account for the 2-body dynamics. Several years ago, Lepage [27 provided a general approach to developing 
effective interactions starting with the now-canonical series of contact interactions. However, in his treatment he let 
these morph into a similar series smeared in configuration space by what amounts to a form factor. The possibility of 
using separable potentials is also mentioned in passing in several of the important papers in the recent history of ET's 
addressing three-body physics including one by Afnan and Phillips |28j . In any case, there is much current interest in 
how best to treat finite range effects that may enter in Efimov physics which is now accessible experimentally in cold 
atom systems where one exploits Feshbach resonances to tune the 2-body scattering length, a, to essentially any desired 
value. As the scattering length diverges, finite range effects become negligible and the physics is universal [29] I30j. 
However, when away from this limit and when interested in how such a limit is approached with a finite, it would 
seem prudent to use an approach that does not implicitly assume range effects are absent at leading order. 

Building on earlier work by Weinberg [31], E. Harms |13j developed a separable expansion for the 2-body potential 
which admits an algebraic matrix solution for the T-matrix. Several authors have studied the convergence of the 
separable expansion of the 2-body potential in 3-body applications 19, 32, 33 , 34 . The separable expansion truncated 
at leading order is historically referred to as the "Unitary Pole Approximation" (UPA). In the present work we focus 
on the ^S'l NN channel and implement the UPA using a dipole form factor with parameters that give a pole in 
the 2-body T-matrix at the 2-body binding energy (of the deuteron) and show how its associated form factor leads 
directly to the 2-body bound state wave function. To establish this methodology in a simple context, we treat the 
2-body channel without tensor coupling and then focus on the nd 5*3/2 (quartet) 3-body channel which allows us to 
avoid complications due to Efimov physics. (Extensions to tensor coupling and the nd doublet case will be treated 
elsewhere.) We then develop a systematic expansion of the 3-body scattering amplitude in inverse powers of the dipole 
range parameter and demonstrate how the ET-ERE based on contact interactions follows as a limiting case of our 
model which further allows us to study its convergence properties in the 3-body context. This exercise demonstrates 
some significant limitations of the ERE in accounting for finite range effects. 

The paper is organized as follows. In Section II we review Harm's separable expansion. We focus on the rank-1 
approximation, the UPA, for the NN problem in the ^Si channel and constrain its single inverse range parameter, (3, 
using the Nijmegen II [3S] phase shifts. In Section III we use the UPA to fix the 2-body input to the 3-body Faddeev 
equations for nd scattering in the ^S3/2 and ^Pj channels. To provide a basis for examining finite range effects, we 
develop an expansion of the resulting Faddeev equations essentially in powers of 7//? where 7 is related to the 2-body 
binding energy, B, via 7 = y/rnB. In Section IV we outline the treatment of the same problem using ET methods. 
We review various common approaches for including effects of the finite range of the 2-body interactions based on the 
Effective Range Expansion (ERE) of p cot (5 for 2-body elastic scattering. We refer to such calculations as ET-ERE 
results. Section V presents comparisons of ET-ERE results at Leading Order (LO), Next-to-Leading Order (NLO) 
and Next-to-Next-to- leading Order (NNLO) with the earlier UPA results similarly expanded in powers of the range 
of the form factor. We find, not surprisingly, that the two approaches are identical at lowest order. Some differences 
appear at NLO but there are qualitative differences at NNLO. Moreover, for s- and p-wave quartet nd scattering, 
we find that both expansions of the 3-body equations suffer some catastrophic failures at NLO and NNLO while the 
full UPA provides a excellent description of the phase shifts in this channel in spite of its simplicity. We present our 
conclusions and future directions in the last section. 



II. SEPARABLE DIPOLE 2-BODY INTERACTION 



As noted in the introduction, separable potentials have a long history in nuclear physics going back to the period 
immediately following Faddeev's breakthrough work on the 3-body problem when, for computational convenience, 
calculations of nd scattering relied almost exclusively on separable NN interactions. The essential feature of these 
approaches that we will exploit here is the explicit inclusion of the range of the 2-body interaction and its close 
association with the bound state wavefunction. Following Harms [13], consider the 2-body Schroedinger equation {h, 
c — 1, ni^^iiifQ'fi — nipj^oton — Tn^ . 



(2) 
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where Gq{E) is the free 2-body propagator. Assume there exists a single bound state a,t E ^ —B. The bound state 
satisfies 

G^\~B)\xb)^V\xb). (3) 
Defining the form factors via the operation of G^^{—B) as 

\i^^)=G-,\-B)\xn), (4) 
we see that the bound state Schroedinger equation can be written as 

FGo(-S)I^b) -AbIVb), (5) 
where A_b = 1. Solutions of the generalization of Eq. [s] 

FGo(-B)|V'«) = A„|V'n), (6) 

form a complete set of kets with normalization condition 

(V„|Go(-B)|7/.„) = -^„,„. (7) 

We call the |'0„)'s the "form factors" of the potential V . (Note that our A„'s are the inverses of those defined by 
Harms.) Also note that, with this normalization convention, the bound state wave function is not normalized. We 
may define the wavefunction renormalization Z via 



d 



Zb' = {xb\xb) = {MG^{-Bf\i'B) = -^^ 



{ijB\Go{E)\i;B}. (8) 

E=-B 



When there is a single bound state, the eigenvalues satisfy 

A„^B < As = 1. (9) 

We observe that the potential can be written as 

V=-T. (10) 

n 

which is verified by substitution into Eq. [6] and then using the normalization condition Eq. [7j After ordering the 
terms in this summation by their magnitudes of eigenvalues beginning with Ai = 1, a "rank-A'^" approximation to the 
2-body potential results from truncating the sum at N terms. The 2-body T-matrix satisfies the Lippmann-Schwinger 
(LS) equation, 

r(p,p-;E^).V(p,p-)+/^ '^'7ry^^ (11) 

J (zn)-^ E+ — 

It is easy to show that 

N 

T(p,p';£;+) = - ^ (pIV-™) A(i?+)„,™ (V'^Ip'), (12) 

n,m— 1 

where 

[A~'(£^+)]n,™ = Xn^5n,m + {-^ n\G q{E+ ^) . (13) 

Thus, we see that the 2-body T-matrix associated with a separable potential to a given rank involves inverting a matrix 
with dimension of that rank. The Unitary Pole Approximation (UPA) corresponds to truncating the expansion of V 
and T at rank-one, including only the form factor for the bound state: 

T/->-|^b) (V-sl, (14) 

and 

T(p,p';ii;+) -> -(pI^b) A(i?+) (V^bIp'), (15) 
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where we have used = 1 and where 



A-'{E+) = 1 + (V's|Go(£;+)|Vb). 



(16) 



To proceed, we make a specific choice for the functional form of the bound state form factor. We select a dipole (or 
Yamaeuchi |36 ) form, namely 



(p|Vs>-AAl/2g(p2) 



(l+p2//3')' 

where g is the form factor which satisfies 17(0) — 1, p'^ — p^, and JV is found by normalizing according to Eq. [t] 



(27r)3 S + p2/„^ 87r(l + 7//3)2 



where 7 = \/mB has been used. Similarly, 

A-\E+) = 1 - 



/3 + 7 



13 + V-mE+ , 

Finally, using the dipole form factor in Eq. [8] gives the wave function normalization factor, 

^ 7(/3 + 7)" 



(17) 



(18) 



(19) 



(20) 



It is clear that, because of the normalization condition, the 2-body T-matrix will have a pole at E — ~B as required. 
It is interesting to note that this choice of form factor means we are, in effect, using the Hulthen parameterization of 
the deuteron wave function [37]. This can be confirmed by using the inverse of Eq. |4| namely. 



{p\xb) = (p|Go(-i?)|Vs), 



(21) 



to find the momentum space bound state wavefunction, {p\xn)- Fourier transforming to configuration space yields 
the radial wave function, 



u{r) — r XB^r) oc (e — e ^^), 
which in turn implies a local r-space NN potential, 

/32_72 



V{r) 



(22) 



(23) 



We are reminded that, in the lore of the day when the Hulthen wave function was in widespread use, the rule of 
thumb was that /? ~ 67 produced a good account of the structure of the deuteron. 



Using Eqs. 15 Q6l and 17 we may construct a fully off-shell elastic 2-body amplitude, /2, via 



hip,p';E)^-^np,p';E)^-^ 



-Mgip") g{p")A{E) 
We then obtain the on-shell elastic scattering amplitude, f(p), and use 

f(p) = f2{p,p;E = p^/m) 

to determine 



1 



p cot S — ip' 



p cot i5 = 



^7(2/3 + 7) 3/32 2/37 -f 72 



1 



2(/3 + 7)2 2/3(/3 + 7)2 



-P 



2/3(/3 + 7)2 



Comparison with the Effective Range Expansion (ERE) 

p cot 6 -■ 



1 '^O 2 , 4 
- + —P +ViP . . . , 

a 2 



(24) 



(25) 



(26) 



(27) 
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allows us to identify 



2(1+7//?)^ _ 3 + 27//? + 7V/32 1 

7(2 + 7//3)' /3(l + 7//3)^ " 2/?3(l + ^/;3)2 ' ^^S) 



Note that the ERE for the dipole form factor terminates at O(p^). To constrain the parameters, we fit our model to the 
phase shifts and binding energy as determined by the Nijmegen II (NIJMII) potential ^SQ which gives 7 — 0.231607 
fm"^ Fitting to the NIJMII ^Si phase shifts for pl„^ < 0.85 fm'^^ we determine /3 = 1.3940 fm"^ Note that 
f3/"/ — 6.019 in keeping with the expectation for the Hulthen wave function mentioned above. We then predict 

a = 5.421 (5.403) fm, tq = 1.772 (1.749) fm and V4 = 0.1357 (0.163) fm^ (29) 

which are in good agreement with; e.g., Babenko and Petrov [35] whose values appear above in the parentheses. 
This level of agreement is gratifying considering the simplicity of the approach. (We observe that a similar level of 
agreement is found upon assuming a Gaussian form factor although the ERE does not terminate at finite order in this 
case.) Note also that the value of (3 is roughly two pion masses which is what we would expect qualitatively for this 



"cutofF parameter; i.e., (5 is "natural" [T]. Finally we observe that none of the quantities in Eq. 28 is ill-behaved as 
13 — > 00; i.e., as the potential becomes a simple contact (i5-function) interaction. In this limit, a — > I/7, ~ (3^^ — » 
and W4 ~ (3~^ ^ as expected. We also emphasize that the ERE is fixed by the on-shell elastic amplitude, Eq. 25 



whereas, as we shall shortly see, the fully off-shell amplitude of Eq. 24 is - in principle - required in the 3-body 
Faddeev equations. 

III. 3-BODY FADDEEV EQUATIONS USING THE UNITARY POLE APPROXIMATION 

Having explicit forms for the off-shell 2-body T-matrices allows one to calculate the 3-body scattering amplitude 
easily. The use of separable potentials in 3-body applications has a long and storied history to which we cannot do 
full justice here. Some recent work which has made use of separable potentials in 3-body calculations can be found in 
Refs. [13 ESI m]. Following Waston and Nuttall [H] but using our sign conventions, the nd ^£"3/2 T-matrix for the 
i*'^ partial wave is given by 

X,{p,p';E)^2Z,{p,p';E)~ — q^dq2Z,{p,q;E) ^{E - ^) X,{q,p';E), (30) 



where A(i?) is given by Eq. 19 identified with the deuteron (or dimer) propagator, and Zi is the particle-exchange 
("ping-pong") amplitude depicted diagrammatically in Fig.[l] 



rn ^^ V>b(|p + p72|)V>b(|p' + p/ 2| 

2 J_i p^ + p'^ — mE+ + pp'x 



ZeXp,p;E)^~Ais- I dx , , — ZZl^^^TZZJZ ^^a^), (31) 



where Pi{x) is a Legendre polynomial and A/5 is the isospin-spin structure factor which is unity for spinless bosons. 
For the nd quartet case, 

A« = C/(i.i.^,l;M)xa(l,l,l,l;0.o).(l)x(-l). ,32) 

where the C/'s are the unitary Racah recoupling coefficients. Note that the form factors for the 2-body bound state, 
ij^B, appear in Z. 

We can define the (fully off-shell) 3-body elastic scattering amplitude for the partial wave by including the 
wavefunction normalization factor, Eq. |20| as 

Tfl 

Mp,p';E) = - — ZbX,{p,p';E). (33) 
ovr 

Then the on-shell 3-body partial wave scattering amplitude is in turn given by, 

Mp) = fi {p, P;E^ (3pV4 - l^)/m.) . (34) 
Using Eq. 30 the Faddeev equation for fi becomes: 

roc 

f,{p,p';E)^f^''"\p,p';E)+ dq q^K,{p,q;E)U{p,p' -E). (35) 

Jo 
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For our form factor model we have 



8A7S7 



(1 + 7//3)' zj,'\p,p';E), 



with the reduced Z, 



(0) 



defined by 



1 



pqx) g{p'^ 

— mE - 



f gV4 - 
pqx 



pqx) 



Pe{x). 



(36) 



(37) 



For e : 



specifically, using Eq. 17 for g{p), we have 



zi,°\p,q;E) 



3pq{p'^ 



g2) 

-1 



In 



p^ + q^ — mE + pq + + <Z^/4 — pq 



+ mi? — 3(7^/4 + — mi? — pq 

+ — mi? + pq 



1 



/32 + mi;-3p2/4 



In 



/32 



74- 



-pq 



q^ — mE — pq /J^ + pp- jA: + q^ + pq 



(38) 



The kernel, Kc, is defined as 



K,{p,q,E) = - 



1 

4A/5 (1 



2Z,(p,g,i?).A(i?-^) 



IT 



zf\p.q;E) 



/3 + 7 



/? + V^3g74 - m£; 



(39) 



It is important to note that the kernel contains the three dynamical elements of the fully off-shell elastic scattering 



amplitude appearing in Eq. 
propagator, A(i? — iq^ /Am) 



24 



namely the two form factors which reside in zf^ as given in Eq 



37 



and the dimer 



The arguments of these functions do not in general satisfy the on-sheUcondition. This 
dependence on the full off-shell structure of the 2-body amplitude is a well known characteristic of the 3-body Faddeev 
equations. As we shall soon see, in the standard ET approach, the form factors are - in effect - set equal to one and 
only the on-shell 2-body amplitude, Eq. |25| enters. 



Using Eqs.[36|and [39| we solve Eq. [35] directly to find the 3-body elastic scattering amplitudes and phase shifts for 
the model. (We employ the method of Hetherington and Schick [39] to avoid ill-behaved integrands in the homogeneous 
term for energies above threshold.) Since 1//3 is roughly the range of the interaction, we can study the effect of the 
finite range explicitly. For this purpose we define Leading Order (FF-LO), Next-to-Leading Order (FF-NLO) and 



35 



m powers 



Next-to-Next-to-Leading Order (FF-NNLO) expressions by expanding f^"™ and Kg on the rhs of Eq. 
of l/f3. The LO expression corresponds to retaining only terms independent of 1//3; this obviously coincides with 
taking the /3 — > oo limit which corresponds to the simplest 2-body contact interaction. In particular the /3 — > oo limit 



of Z, 



(0) 



gives. 



Z^°\p,q,E) 



1 



In 



P 



q^ — mE + pq 



(40) 



pq \p^ + q"^ — mE — pq ^ 

a factor that is ubiquitous in ET-ERE formulations. It is gratifying that this limit can be taken straightforwardly 
without relying on any particular prescription for renormalization of divergent integrals. We note that the /? — *■ oo 
form of Zq°' is the only form which appears in the ET-ERE 3-body equations at any orderl (See, for example, 
Ref. [H], specifically the discussion following their Eq. 34.) In contrast, we find, upon expanding z'^\p,q,E), 



zl^\p,q;E) 



1 



■In 



p^ 



mE 



pq 



mE 



pq 



1 



3(p2 + g2) _ 



4/32 



1 

4/32 



(41) 



00 form starting at NNLO. We will comment on the "extra" 



pq \p^ + q^ 

which implies this quantity is modified from its f3 
terms in the discussion to follow. 

We next give explicit formulas for the s-wave amplitude; the subscript £ = will be implicit henceforth. Expanding 
the s-wave Born term, Eq. |36| to second order, using Eq. |38| we find. 



fE 



8A/S7 f 1 i^fP^ + (1^ ~ + PI 
3 \ pq + q^ — mE — pq 



[pq 

/3 4/32 



(3p^ + 3g2 + 1272 - 8mE) 




(42) 
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In similar fashion we expand the s-wave kernel, Eq. 39 to find, 



Kip,q;E) 



1 



f 1 i^f 1^ ^ ^"-^ + PI 



TT ^3(/2/4 - mE [ pq VP^ + - rriE - pq 



i + ^h + V^qy^-mE) 

+ (I2p2 + I5q2 + _ 3g^£; _^ 407v/3g2/4 - mE) 



13^ 



(43) 



Similar expressions can be obtained for the higher partial waves. 
We immediately see that the constant term, —4/(3"^, in Eqs. 



42 



and 



43 will give a contribution at NNLO which 
will have no counterpart in the ET-ERE NNLO expressions (see Eq. 41). It is interesting to note that the form 
of this term is identical to the lowest order 3-body contact interaction invoked in the ET-ERE formulation. Upon 
reflection, the appearance of such a term is not surprising. We know that, for the 3-boson case (A/g = 1), the 3-body 
amplitude shows a log-periodic dependence on the upper limit of the q integration in the homogeneous term |10l 111] . 
Regularization is accomplished by introducing a simple 3-body contact interaction like the — 4//?^ term just mentioned. 
In the UPA version of the 3-boson problem, the form factor regulates the integral so no log-periodic dependence on 
the integration limit is found. Therefore, when ET-ERE formulations add 3-body terms, at least in part, they are 
compensating for additional 2-body finite range effects not accounted for in that approach. 



Combining Eqs. 42 and |43[ we find the Faddeev equations for our model at the various orders in the range of the 
interaction, 1//3. First, at leading order, we find 



fLo{p,p';E) = 



8Ajs7 ^ J_ 
3 pp' 



In 



p- 



q^dq — In 
pq 



p — mE 

2 



PP_ 

— mE — pp' 



p^ -\- q^ — mE + pq 



mE - 



pq 



fLo{q.p';E) 



-7+ ^iq^l-i - mE' 



(44) 



which is, not surprisingly, identical to the ET-ERE equation at LO which will appear below (hence no notational 
distinction between models is needed at LO). We again emphasize that this expression corresponds to the limit 
/? — > oo which implies a simple 2-body contact interaction and that passing to this limit is trivial in the present 
formulation. With A/5 1, this expression corresponds to the Skornyakov-Ter-Martirosyan (STM) equation for three 
identical spinless bosons [1CT1I1T| . and identical to the ET-ERE LO Faddeev equation (TUl HUIISI |31] (and many more; 
see,e.5., Refs.[29, 3^ for extensive reviews of this subject). 

Similarly, the NLO Faddeev equation for the form factor model is readily determined to be 



f!,loM,E) = 



3 

, 2A/S 



1 

pp' 



In 



p^ 



pZ _|_ pi 

1 



p — mE 
2 



q^dq — In 
pq \P' 

f^lo{q,p'-E) 
-7+ ^iq^/A-mE' 



1P1P_ 

— mE — pp' 



p^ + q^ — mE 



-pq 



mE 



pq 



1 + ^(7+ mi?) 



(45) 
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Likewise for the FF-NNLO Faddeev equation, we find 



8A/S7 I 1 fp'^ + p'^ - mE + pp' 



PP 



p^ + p' — mE — pp' 



1 + 3^ + 31, 



1 

4/3- 



+ 



) 

q^dq 



In 



pq 



mE + pq 



mE — pq 



1 + + \/W/4^^mE) +\7T7sh+ V3gV4 - mE) 



2/3 



~ ~ ^^^^ ^ V3gV4 - mE )) 

„ fNNLoil^P''^E) 

-7+ ^/3g2/4 - mi;' 



4 



(46) 



where the first three terms in square brackets in the kernel are organized to resemble the corresponding expression 
for the ET-ERE to be considered next. 



IV. 3-BODY ET-ERE FADDEEV EQUATIONS 

While the finite range expansion for our form factor model is straightforward, extending ET-ERE formulations to 
include corrections due to the finite range of the 2-body interaction in 3-body applications is not trivial and has been 
the subject of considerable attention in the literature; see Refs. [HllIS] for extensive discussions of this treatment and 
for additional related references. 

To emphasize the relation between our FF calculations as formulated in the preceding Section and the ET-ERE 
treatment, we obtain the latter by beginning with Eq. [30] our full FF 3-body Faddeev equation for the amplitude X 
which we re-state as: 

Xe{p,p';E) = 2Zeip,p';E)^ ^ q^dq 2Zeip,q;E) a(^E - X,{q,p'-E), 



where A(i?) is is given by Eq. 16 and can be identified with the deuteron (or dimer) propagator and where, for s-wave 



scattering and assuming is real, is given by Eq. |31|as 



Zq{p,p';E) 



-Ajs- J_^dx 



A m 

= -hjs—M 



2 

m 



= -his'^n Z'^'>{p,p';E), 



' ^ j{\p+p'mg{\p'+pm 

I p'^ + p'^ — mE + pp' X 

r(0) 



(47) 



where Af is the normalization factor given in Eq. 18 and g is the form factor (see Eq. 17) 



We proceed toward the corresponding ET-ERE expression by first examining the homogenous term which we can 
write as 



mAjsAf 
27r2 



q dq 



p. 



q^ — mE + pqx 



3q2 
4to 



)-X{q,p'-E), 



(48) 



where the ^ = subscript has been suppressed. The next step is not physically justified but amounts to a prescription 
which yields the ET-ERE expressions. The arguments of the form factors are set to mE — iq^ /A in which case they 
no longer depend on the angular variable x. Thus the angular integral can be performed in the usual way and the 
homogeneous term becomes 



2A 



IS 



q^dq 



1 

pq 



In 



p- 



mE + pq 



mE ■ 



pq 



mAf 
in 



g{mE-3qyAf -AiE-^) 



X{q,p'-E). 



(49) 



9 



Now, according to Eqs. [T2| 13 and 24 the qu antity in the square brackets is precisely the on-shell 2-body elastic 

^ mE — 3(7^/4. Now using, Eq. 



scattering amplitude, f(p), evaluated at p 
Eq. [27) we can rewrite the homogeneous term as 



25 



and the Effective Range Expansion, 



c^dq — In 
pq 

q^dq — In 
pq 



p- 



pq 



mE 



pq 
pq 



r\p= ^mE-3qy4) ■X{q,p'-E) 

X{q,p';E) 



p^ + q^ -mE-pqJ ^1+ 'f{mE~3q'^/A) + V4,{mE~3q'^/A)^--- + ^iq^j^-mE 



(50) 



Next we look at the Born term in the equation for X. Using the expression for Ze{p,p';E) appearing in Eq. [31] 
the Born term in Eq. [30] for our full FF 3-body Faddeev equation for £ — is, 



2Z{p,p';E) = -mAis J dx 



V's(|p + p72|)^i3(|p' + p/2|) 



p^ + p 



i'2 



iE+ 



pp'x 



where here and henceforth the I — Q subscript is suppressed. To make contact with the ET-ERE approach, we ignore 
the momentum dependence of the V's's, perform the integral over x and rewrite the Born term as 



2Z{p,p'-E) 



-m Ajs ■Af ■ — 7 In 
PP 



7 1 /• 

p +p 



mE 



PP 



pi _ jj-i^JTj _ ppl 



pz 



(51) 



where J\f is the form factor normalization (see Eq. 18 1. We remind the reader that the assumption made above in 
deriving this simplified Born term ~ ignoring the momentum dependence at the vertices of the "ping-pong" diagram 
- is always made implicitly and its validity is never checked. We next find the Born contribution to the 3-body elastic 
scattering amplitude, /, by making use of the relation between / and X as given in Eq. [33] 



ET-ERE 
Born 



m^Ais 



tERE 



37r 



1 

PP' 



In 



p^ + p' — mE + pp' 



(52) 



where Z^^^ ^ ZbM is the ET-ERE wave function renormalization to be discussed shortly. 

Then, using the expression in Eq. 52 for the Born term and the expression for the homogeneous term in Eq. 50 the 
ET-ERE Faddeev equation for the 3-body scattering amplitude becomes 



fET-ERE 



7ERE 



37r 
2 

+ - 



1 

pp' 
1 



In 



P 



q^dq — In 
pq 



pi _|_ pi 
p^ + q 



p — mE 

2 



PP_ 

— mE — pp' 



mE - 



pq 



q^ — mE — pq ^ 

jET-ERE(^q^pi.E) 



+ ^{mE - 3(?2 /4) + Vi{mE - 3q'^ /A)^ ■ ■ ■ + ^3(?2 /4 - mE ' 



(53) 



It turns out in practice to be more convenient to use a form of f ^ (p) which is equivalent to using a v ersion of the 
ERE based on an expansion of pcot 6 about p^ = —7^ rather than about = as is implicit in Eq. 27 Thus 



r\p)--7+^(p^ + 7^) + -;(/ + 7Y 



ip, 



(54) 



where there is clearly a simple relation between the original ERE parameters in Eq. [27] and their primed counterparts 
in the last equation. (For example, we have tq = 1.772 fm but Tq = 1.747 fm based on our dipole form f acto r.) We 
will adopt this latter prescription to generate our ET-ERE NLO and N^LO Faddeev equations from Eq. [44] Other 
methods such as the ^-expansion approach appear in the literature; we refer the reader to Refs. |12j for details. In 
all ET-ERE cases we have examined, the various forms of the ERE expansion yield results which are qualitatively 
similar. 
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Using the form of f(p) given in Eq. 54 our putative higher order ET-ERE s-wave Faddeev equation is then 



OTT 





^ERE ^ 


In 




vv' 




^2A,s 


poo 1 


In 


TT 


Jo PQ 





p2 _|_ _ ^^^^ _|_ 
p2 _|_ p/2 _ ^^^^ _ 



r' 



• • • + ^3974 - toS 



(55) 



which apparently includes contributions from the ERE to all orders. However, as has been noted frequently in 
the literature [TTJ [23], use of this form of what amounts to the dimer propagator in 3-body Faddeev equations is 
problematic as it introduces spurious poles in the integrand of the homogeneous term which lead to unphysical effects, 
the most serious of which are imaginary phase shifts that are non-zero below the breakup threshold or that are 
negative. A nearly universally applied prescription for avoiding (some of) these pitfalls is to use (see, e.g., Ref. |44j 
for one of the earliest discussions of this point) for f{p), instead of Eq. [54] 



f(p) 



-7 — ip 



N„ 



E 

n=l 



(56) 



which obviously only takes into account the 0{p^) term in the ERE but which preserves the pole structure of the LO 
integrand (Eq. 44 1 in the homogeneous term. We henceforth adopt the terminology frequently used in the literature 
to designate calculations employing Eq. 56 as NLO or N^LO for Nmax = 1 or 2, respectively. We adopt these 
prescriptions to generate our ET-ERE NLO and N^LO Faddeev equations by using the form of f(p) appearing in 
Eq. [56] in the 3-body Faddeev equation shown in Eq. [55j 

To find explicit ET-ERE 3-body equations, it remains to specify the wave function rcnormalization, Zere- Fol- 
lowing, e.g., Ref. [I0j[nj|42] , we write 



'ERE 



m d 
'i^ dE 



iEREiP = ^^E+). 



(57) 



E=-B 



According to Eq. 



56 



at LO, islAP = ^rnE+) 



—7 + \/—mE+ which implies, with B = 



7ERE 



8717 

nrrt 2 



(58) 



Similarly, we readily find 



yERE 
^NLO 



yERE 
^NNLO 



^(i + ^W + ^M- 



(59) 



Thus, at NLO we find 



lET-ERE ( I 77i\ 
Inlo {P^P\E) 



87A/5(l + r(,7) 1 ^ (p^+p 



In 



2A 



IS 



PP' 

q^dq — In 
pq 



mE - 



PP 



p^ + p' — mE — pp' 

2 , ^2 I „„\ f ET-ERE 



iq,p';E) 



q^ - mE - pq) + ^3g2/4 - mE 



+ ^{l+^JWI^-mE) 



(60) 
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while at N^LO 



Jnnlo iP'P'^) 



3 pp' 

„2 , 



p'^ — mE + pp' 



2A 



IS 



^p^ + p' — mE — pp' ^ 
'■^-mE + VQ\ f^ll^''''iq,p';E) 



q^dq — In ( ^ , 

pq \p'' + q^ ~ mE - pqj + y'Sq'^/A - mE 



+ (j{l+V5qy^-mE)^ 



(61) 



Comparing the FF and ET-ERE Faddeev equations at NLO, we find Eqs . [45] and [60| would be identical if we could 
equate Tq with 3//3. This relation is found to be satisfied numerically with reasonable accuracy (t'q ~ 1.75 fm vs 
3//3 ~ 2.16 fm), but the difference is a real one since we must use the parameters as set by the 2-body phenomenology 
if we are to be consistent. We then must conclude that incorporating finite range effects via our separable potential 
with a simple dipole form factor gives rise to contributions to 3-body scattering which cannot be extracted exclusively 
from the 2-body ERE. 

Next we compare the FF and ET-ERE Faddeev equations at NNLO; i. e. Eqs. 46 and 61 where significant differences 
are apparent. The three terms in the square brackets in the homogeneous term of the ET-ERE Faddeev equation, 
Eq. [61] also appear in that position in the FF NNLO equation, Eq. ]46] These terms can be traced back to the 
expansion employed in the ET-ERE for f(p) which is specified in Eq. 56 However, extra terms at 0{1/ 0^) are present 
in the FF expression including the term which resembles a 3-body contact interaction as discussed previously 

(see Eq. 41 1. We conclude that, while the ET-ERE, as implemented employing Eq. 56 indeed avoids pathological 



pole structure of f(p) in the homogeneous term of the ET-ERE 3-body Faddeev equations, it ignores the form factors 
appearing in Fig. [l] and hence fails to account for the momentum dependence contained in the form factor vertices 
associated with the 2-body interaction. Consequences of these shortcomings will be examined in the next Section. 



V. RESULTS 



We now present calculations based on the expressions developed thus far. As discussed in Sect. ]lT] our 2-body 
input is constrained by the Nijmegen II (NIJMII) np phase shifts in the channel. We first determine pcoi8 and 
7 = 0.231607 fm^i from the NIJMII potential. We then adjust /3 in Eq. [26]to give a fit to pcotb up to ^ 0.85 
fm^2, finding — 1.3940 fm~^. We then use the expressions in the previous section to determine s-wave phase shifts 
for the ^53/2 (quartet) channel in nd elastic scattering. To provide a context for comparison, we also show 3-body 
phase shifts from full Faddeev calculations using realistic NN potentials, Refs. [151 H5] . 

For our first set of calculations we demonstrate the convergence of the finite range expansion by artificially letting 
/3 ^ 2 X /3 ~ 2.79 fm~^ which means the 2-body effective range parameter, rg, decreases by roughly one half. The 
expansions of the 3-body Faddeev equations which give our Leading Order (LO), Next-to-Leading Order (NLO) and 
Next-to-Next-to-Leading Order (NNLO) expressions of the preceding Section - which are expansions in 7//? in the 
FF case or '^t'q for the ET-ERE case - will converge rapidly since, eg, 7//? ~ 1/6 ^ 1/12. Our calculations of the real 
and imaginary phase shifts for this artificial case appear in Figs. [2] and [3] respectively. 

For the real part of the phase shifts. Fig. J2] both the FF and ET-ERE calculations have converged at NLO and 
that there is little difference between the two approaches here. (Agreement with the full Faddeev calculations is 
understandably poor since we have artificially increased /3.) Examination of the imaginary part of the phase shifts. 
Fig. ]3] is a bit more revealing. Here we see some differences at the various orders of the expansion, with the NLO 
results lying below the NNLO and full FF values. (The fact that the LO calculations are in near agreement with the 
NNLO and full FF results is coincidental.) At each order (LO, NLO and NNLO), the FF and ET-ERE calculations 
are quite similar suggesting that, in this artificial case where convergence has been accelerated, the two approaches 
are nearly equivalent. Closer inspection reveals two interesting features. First, both NLO calculations yield slightly 
negative imaginary phase shifts from threshold up to about p = 110 MeV/c. Though these negative imaginary phase 
shifts are small in magnitude, the are large compared with our numerical uncertainties and, from a physical point of 
view, they represent a serious failure as they are incompatible with unitarity. This failure is evidently corrected in 
this case by going to NNLO where the imaginary phase shifts are positive above threshold as we would expect. The 
second point of interest is that, at NNLO, the FF phase shifts arc much closer to the full FF values than are the 
ET-ERE results. This is perhaps not surprising since the FF NNLO is based directly on the expansion the full FF in 
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powers of 1//? and presumably includes, e.g., some of the off-shell properties of the 2-body T-matrix which cannot be 
accounted for in the ET-ERE formulation. 

We now turn to realistic nd calculations which employ the value oi (3 — 1.3904 fm^^ fit to the NIJMII phase shifts. 
We first consider the quartet scattering length ^03, values for which are presented in Table |l] Here our results for the 
various calculations under discussion are compared with the old experimental results of Dilg et al. |47j and the modern 
full Faddeev results of Witala et al. [48 . We focus on the latter result as it corresponds to the specific case of the 
NIJMII potential without electromagnetic interactions and hence most closely resembles our NN input. We see that 
the LO, NLO and NNLO calculations of both types steadily converge toward the Witala et al. value of 6.325 fm, but 
that the FF results appear to converge somewhat faster. Our full FF calculation agrees almost exactly with Witala 
et al. value. Next we look at the real and imaginary s-wave phase shifts which appear in Figs. |4] and |5] respectively. 
For this realistic case, differences among the various calculations are apparent for both the real and imaginary phase 
shifts. For the real phase shifts, both the FF and ET-ERE calculations appear to have converged at NNLO, but the 
convergence of the former seems a bit faster. Both the FF NNLO and full FF calculations agree with the full Faddeev 
results quite well while the ET-ERE results are arguably just as good. 

Examination of the imaginary phase shifts in Fig. [5] leads to rather different conclusions. First, the failure in 
the form of negative imaginary phase shifts has become catastrophic. The NLO and NNLO calculations of both 
types yield large negative imaginary phase shifts over the entire momentum range covered. We have carried out 
higher order calculations (not shown) and find that the failure is corrected by N^LO and N'^LO. Moreover, the FF 
NNLO calculation again seems to be appreciably closer to the full calculations than its ET-ERE counterpart. These 
findings suggest (i) that finite range expansions of either type should not be relied upon without carefully checking 
convergence properties and (ii) that expansions based on the on-shell information contained in the 2-body Effective 
Range Expansion may be inadequate considering the appearance of the full off-shell T-matrix in the original Faddeev 
formulation and in the FF treatment considered here. We also note that the full FF phase shifts are in excellent 
agreement with the full Faddeev calculations shown. 

As an additional test of the finite range expansions, Table [ll| shows calculations of the (center of mass) total, elastic, 
break-up, and back-angle cross sections for quartet nd scattering at Eiab — 10 MeV, summing partial waves up to 
I — ^ for each of the models under study. For comparison purposes the table also shows the calculations of Koike 
and Taniguchi [TS]. (These authors point out that the nd elastic cross section at 180° is dominated by the quartet 
contribution so separation from the doublet is not required.) The ET-ERE and FF NNLO calculations are typically 
close to one another, but each falls appreciably below both the full FF calculation and the realistic calculations of 
Koike and Taniguchi. The discrepancy of the NNLO calculations with the other results for the breakup cross section 
is traceable in part to the negative imaginary phase shifts of the former. It is gratifying to note that results from our 
simple FF model are consistently within 10% of the Koike and Taniguchi values. 

We now turn to calculations of the i = 1 nd quartet phase shifts shown in Figures [6] and [7] and similarly show LO, 
NLO, and NNLO calculations for ET-ERE while taking the FF to N^LO as well as the full form factor calculation 
and the full Faddeev calculations of [45] . 

For the real part of p-wave quartet phase shift. Fig. |6] we see that (i) the full form factor calculation is in good 
agreement with the lower range of the full Faddeev calculations, (ii) the ET-ERE NLO and FF NLO are essentially 
identical but differ at NNLO. We have carried out higher order calculations (not shown) which oscillate about, but 
converge to, albeit slowly, the full FF result. 

The situation is more dramatic for the imaginary part of the p-wave phase shift as shown in Fig. [7j Again the 
full form factor calculation is in fair agreement with the full Faddeev calculations. The LO calculations look quite 
good as do both the ET-ERE NLO and FF NLO calculations which are essentially identical. However, the ET-ERE 
NNLO and FF NNLO diverge from one another, apparently due to the extra terms in the latter (Eq. [46]) which are 



not present in the former (ET-NNLO expression, Eq. 61 1. Again, the higher order FF calculations oscillate about 



the final full form factor results. Generally, the calculations relying on an expansion in the range of the interaction 
display poor convergence properties while the full form factor results are well behaved and, indeed, give an excellent 
account of both the real and imaginary phase shifts examined here. 

At this point, we address two sets of ET-ERE calculations which do not use the prescription implied by Eq. |56[ 
The first is the work of Bedaque and Griesshammer [43 in which they present calculations of quartet s-wave phase 
shifts based on the ET approach. These constitute, to the best of our knowledge, the first ET evaluation of the 
quartet phase shifts above the breakup threshold. Overall, they they find good agreement with the results of full 
Faddeev calculations. Their most ambitious calculations include explicit pion loop contributions which precludes any 
comparison with the present work. They also show what they refer to as "perturbative NLO calculations with pions 
integrated out" . We have not established how to compare these results with the present work. However, they also 
show "pionless NNLO" calculations - based on a formulation appearing in Ref. [9] - which correspond to using our 
Eq. 55 after setting v'^ and all higher order terms to zero. Our attempts to reproduce these specific calculations reveal 



that they are subject to problems arising from spurious poles as discussed in connection with Eqs. 55 and 56 in the 
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preceding Section and as noted in Ref. (5], too. Our comparable s-wave results are very similar to those of Bedaque 
and Griesshammer [43] and, in fact, are quite similar to our FF calculations appearing in Figs. |4] and [s] where we 
show the real and imaginary parts of the phase shifts, respectively. As these figures show, such calculations agree 
well with full Faddeev calculations. However, for this case we calculate small negative imaginary s-wave phase shifts 
below threshold (they are much larger than our numerical uncertainties) in violation of unitarity. Moreover, using the 
same approach to calculate p-wave phase shifts results in catastrophic failure in the form of large negative imaginary 
phase shifts at all energies. Both problems depend critically on the upper limit of the momentum integration in the 
relevant Faddeev equation. Since the spurious pole referred to above appears at k ~1 fm^^, the problem with negative 
imaginary phase shifts disappears if the upper limit of momentum is less than 1 fm"^ so that the "bad pole" cannot 
contribute. This dependence on the upper cutoff on momentum is undesirable and, in any case, is consistent with the 
discussion of the limitations of ET calculations as discussed above in connection with Eqs. [55] and [56] 

The second of these papers which present ET-ERE calculations which do not rely on the prescription implied by 
Eq. |56]is the work of Griesshammer in Ref. [TT] which makes use of the so-called "Z-expansion" method, originally 
suggested in the 2-body context by Phillips, Rupak and Savage [49]. The essential feature of Griesshammer's method 
is to use (see Eq. 2.12 in Ref. [IT), instead of the fully summed form of Eq. 55 



1 



1 



-7 - 1 _ 



{Zb-1) 



(62) 



the expanded form of Eq. 56 



-7 — ip 



N„ 



E 



{Zb - 1) 
27 



(-7 - ip) 



(63) 



where Zb is the wave function renormalization appearing in Eqs. [20] and |59[ The Z-expansion quartet s-wave 
imaginary phase shifts (Griesshammer does not show these phase shifts) at NLO and NNLO using Eq. [63] oscillate 
wildly and are negative over much of the region from just above breakup to p= 200 MeV. Results using the fully 
summed form of Eq.[62]are satisfactory and quite similar to our FF calculation s ap pearing in Figs.[4|and[5]pro?;irferf the 
upper limit of momentum is less than 1 fm^^ so that the spurious pole in Eq. [62[ cannot contribute. Otherwise, small 
negative imaginary phase shifts are again observed below threshold. Our calculations for the Z-expansion quartet 
p-wave phase shifts at LO, NLO and NNLO using Eq. [63] agree well with those appearing in the Griesshammer paper 
and are consistent with our Figs. [6] and [7] However, the fully summed imaginary phase shifts are again large and 
negative at all energies unless the upper limit of momentum is less than 1 fm~^ so as to eliminate the spurious pole 
contribution. 

In closing this Section, we conclude that the ET-ERE approaches we have discussed in this paper, which ignore 
the finite range of the 2-body interaction at the outset and then try to accommodate it approximately via expansions 
based on Eq. [56j face a difficult challenge. We conclude that, in this regard at least, our FF formulation provides a 
superior starting point for the treatment of 3-body physics. We observe that finite range effects are relatively large 
for the case of quartet nd scattering considered here since the expansion parameter ^tq ~ 37//? ~ 0.45 is relatively 
large. In the doublet channel with both singlet and triplet contributions, the same physics occurs in the NN triplet 
contribution, but there is also a comparable contribution from the NN singlet channel where 7ro — > ro/a2 ~ —0.12. 
Even more advantageous is the situation for atomic helium where r^lai ~ 0.075. However, even in these cases, 
finite range effects can be important and hard to account for using, e.5., the ET-ERE prescription when calculating 
observables such as 3-body recombination coefficients which necessarily involve momentum transfers large enough to 
be above the breakup threshold and which depend on the imaginary part of the elastic phase shifts |26j . 

It also important to observe that, in spite of the significant differences between them, both the full FF and the 
ET-ERE calculations use the same on-shell 2-body input, namely parameters fit to the deuteron binding energy and 
low energy phase shifts. However, the three-body kernel probes the 2-body amplitude off-shell. Specifically, the form 
factor treatment provides momentum dependence in the "blobs" of the ping-pong diagram (Fig. [ij which softens the 
integral in a manner consistent with the 2-body phenomenology which results in a more accurate treatment of the 
off-shell behavior of the 2-body T-matrix than does the ET-ERE formulation without such regulating terms. We 
have examined alternative parameterizations for the form factors, such as Gaussians [33] [23] [51] , and found that the 
specific analytic form of the form factors - and the off-shell behavior they imply - is not important provided the 
2-body phenomenology is respected. 
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VI. SUMMARY AND CONCLUSIONS 



In this work we have studied the convergence properties of finite range expansions of 3-body scattering ampHtudes 
by starting with a simple separable form factor model. The full (unexpanded) form factor calculations give an excellent 
account of 3-body nd s- and p-wave phase shifts and other observables, too. (For simplicity, we have ignored 

the channel in the present paper. Calculations to be presented in a future publication show that full FF doublet 

calculations similar to those presented here for quartet scattering are in similarly good agreement with the relevant 
full Faddeev calculations.) We developed a systematic expansion of the 3-body scattering amplitude in the inverse 
of the form factor parameter, /?, which is proportional to the range of the 2-body interaction. The resulting 3-body 
equations were referred to as FF-ERE expressions. We then sketched a derivation of the standard Effective Theory 
form of the 3-body Faddeev equations in which range effects are included using the Effective Range Expansion for 
the 2-body scattering. We were then able to compare directly these two expansions at Leading Order (LO), Next-to- 
Leading Order (NLO) and Next-to-Next-to-Leading Order (NNLO). Perhaps the most significant difference between 
the two approaches is that the ET-ERE formulation lacks the momentum dependence in the "blobs" in the ping-pong 
diagram (Fig. [T]) arising from the form factors. 

Comparisons of numerical results with full Faddeev calculations using realistic modern NN interactions revealed, 
first of all, that our full Form Factor calculations agree very well with the former results in spite of the simplicity 
of our dipole form factor. Both expansions, namely the FF-ERE and the ET-ERE, showed serious deficiencies, the 
most worrisome of which were large negative imaginary phase shifts at NLO and NNLO, a situation incompatible 
with unitarity. This particular shortcoming was only overcome in the FF-ERE case by going to N'^LO and even 
then, agreement with the full Faddeev results was poor. Poor convergence for both expansions was observed for 
p-wave phase shifts, too. We conclude that approaches which do not treat 2-body finite range effects directly and 
employ expansions such as those used in the ET-ERE or FF-ERE may be unreliable. Moreover, the technical problem 
frequently encountered in any attempt to include 2-body range effects to all orders, namely the appearance of spurious 
poles in f(p) - as discussed in detail in Section 3.4 of Ref. [LT\ - never appears in the full FF formulation as it is simple 
to show that there is only one pole corresponding to the single bound state found for any separable 2-body potential 
of a reasonable form. 

The success of our simple full Form Factor results, especially in the regime where the standard ERE expansion 
breaks down, is gratifying and suggests that the approach warrants further study. We first note that, upon writing 
our dipole form factor as a geometric series 

1/(1 + qyp^) = 1 - - _ . . . , (64) 

we see that it automatically sums a certain set of contact terms to all orders! It is therefore perhaps not surprising 
that our simple FF calculations behave so differently from those based on the expansions discussed above. 

In the context of our Form Factor calculations, there is a fundamental question of whether the simple single form 
factor method (the Unitary Pole Approximation) we have employed here can be established as the first step in a 
sequence of controlled approximations. Preliminary work strongly suggests that we can answer this question in the 



affirmative. In a separate paper, we will use the Harms' expansion of an arbitrary potential via Eq. 10 namely, 

N 

V = -J2 l^n>A„(V'„l, (65) 

n=l 

to construct a systematic method of improving upon the UPA by using separable potentials of successively higher 
rank, a procedure facilitated by the fact that the Harms expansion automatically orders terms by their importance at 
low momentum scales. Such a procedure is satisfactory from a fundamental point of view since the Harms expansion 
of the potential given above is exact in the limit that the sum includes all eigenkets of Eq. [6j The practical utility 
of the method depends, of course, on the convergence properties of this expansion. Again, preliminary work shows 
that, e.g., for realistic potentials describing the ^Sq NN channel, only two or three terms of the Harms expansion are 
required for a very accurate reproduction of the on- and off-shell T-matrix for momenta less than ~ 2 fm^^. This 
convergence can be further enhanced by making Harms' expansions of Viow's determined from realistic NN potentials 
using, e.g.. Similarity Renormalization Group methods (SOJ [5TJ [SU [S5]. We will also have to extend our method to 
treat the tensor interactions in the full ^Si-^Di channel. Moreover, our 3-body calculations using these separable 
interactions must be extended to the nd doublet channel and to the case of three identical spinless bosons {e.g., "^He 
atom-dimer scattering) where Efimov physics is present. All of these points will be addressed in a future publication. 
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Model 


"as (fm) 


Experiment 
Full Faddeev 
LO 

ET-ERE-NLO 
ET-ERE-NNLO 
FF-ERE-NLO 
FF-ERE-NNLO 
Full FF 


6.35±0.02 

6.325 
5.091 (19.7) 
5.935 (6.3) 
6.201 (2.2) 
6.089 (4.0) 
6.280 (0.9) 
6.322 (0.05) 



TABLE I: Calculations of nd quartet scattering length, ^as, for various calculations are compared with an (old) experimental 
value [47! and results of full Faddeev calculations using modern NN interactions gS]. Percent differences relative to the full 
Faddeev value appear in parentheses. 



Model 


atot (mb) 


(Tei (mb) 


aBR (mb) 


1^(^ = 180°) (mb) 


ET-ERE-NNLO 
FF-ERE-NNLO 
Full Form Factor 
Koike &Taniguchi 18 


694.8 (85.5) 

684.9 (84.3) 
733.8 (90.3) 

812.3 


673.9 (87.5) 
652.4 (84.7) 
695.2 (90.3) 
770.2 


20.8 (49.4) 

37.5 (89.1) 

38.6 (91.7) 
42.1 


152.3 (82.3) 
160.8 (86.9) 
172.1 (93.0) 
185.0 



TABLE IL Calculations of center of mass total, elastic, break-up, and back-angle cross sections for quartet nd scattering at 
Eiab = 10 MeV, summing partial waves up to £ = 8. Percentages of the Koike and Taniguchi [TH] values are in parentheses. 
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FIG. 1: (Color online). Diagrammatic representation of the Born ("ping-pong") term of 3-body particle-dimer elastic 
scattering Faddeev equation. The " blobs" are the form factors which depend on the magnitudes of the relative 2-body 
momenta which are |p/2 + q] and |p + q/2| for the lower and upper form factors, respectively. 
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FIG. 2: (Color online). The real part of L = nd quartet phase shifts for the artificial case, /3 —> 2/3. The solid curve is the 
full FF calculation, Eqs. |35|39[ and the dotted curve is the LO calculation, Eq. |44| The other curves are NLO calculations 
in the ET-ERE formulation (dash), Eqs. pOpl] and FF-ERE formulation (dash-dot), Eqs. |45|46| The NNLO calculations are 
identical to the NLO indicating rapid convergence as expected for this case using an artificially short range parameter. All FF 
calculations are based on a dipole form factor with with an original inverse range parameter of {3 — 1.394 fm~^ determined by 
fitting to pcotS as given by the Nijmegen-II ^Si phase shifts |35| . In the ET-ERE calculations, multiplying the original /3 by 2 
implies an effective range parameter ro ~ 0.97 fm. Also shown are the full Faddeev calculations of Huber, et al. [45] (triangle) 
and Kievsky, et al. [IS] (plus). 
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The other curves are NLO and NNLO 



FIG. 3: (Color onhne). The imaginary part of L = nd quartet phase shifts for the artificial case, (3 ~* 2(3. The solid curve is 
the full FF calculation, Eqs. |35|39l and the dotted curve is the LO calculation, Eq. 
calculations in the ET-ERE formulati on (das h and dash-dot-dot, respectively), Eqs. 
and dash-dash-dot, respectively), Eqs. 45146 



60|61[ and FF-ERE formulation (dash-dot 
All FF calculations are based on a dipole form factor with with an original inverse 
range parameter of /3 = 1.394 fm~^ determined by fitting to pcotj as given by the Nijmegen-II ^Si phase shifts ^35^. In the 
ET-ERE calculations, multiplying the original /3 by 2 implies an effective range parameter ro ~ 0.97 fm. Also shown are the 
full Faddeev calculations of Huber, et al. [l^ (triangle) and Kievsky, et al. [46] (plus). Note that breakup threshold is at 
p = 53 MeV/c. 
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FIG. 4: (Color online) 

and the dotted curve is the LO calculation 
and dash-dot-dot, respectively), Eqs. |60|6T 



The real part of L = nd quartet phase shifts. The solid curve is the full FF calculation, Eqs. |35|39| 

The other curves are NLO and NNLO calculations in the ET-ERE (dash 



Eq. 



44 



and FF-ERE (dash-dot and dash-dash-dot, respectively), Eqs. 45|46 formulations. 
All FF calculations are based on a dipole form factor with (3 = 1.394 fm"'^ determined by fitting to p cot 5 as given by the 



Nijmegen-II ^Si phase shifts |35j. In the ET-ERE calculations, this corresponds to an effective range parameter 
Also shown are the full Faddeev calculations of Huber, et al. (triangle) and Kievsky, et al. 46 (plus). 



1.747 fm. 
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FIG. 5: (Color online). Same as for Fig. [4] but for the imaginary part of the nd quartet phase shifts. Note that breakup 
threshold is at p = 53 MeV/c. 
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FIG. 6: (Color online). Same as Fig. |4] except for the real part of L = 1 nd quartet phase shifts. 
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FIG. 7: (Color online). Same as for Fig. [4] except for the imaginary part of the L — I nd quartet phase shifts. Note that 
breakup threshold is at p = 53 MeV/c. 



